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Abstract. 

We consider a Jepsen gas of N hard-point particles undergoing free expansion on 
a line, starting from random initial positions of the particles having random initial 
velocities. The particles undergo binary elastic collisions upon contact and move freely 
in-between collisions. After a certain ordering time T Q , the system reaches a "fan" 
state where all the velocities are completely ordered from left to right in an increasing 
fashion and there is no further collision. We compute analytically the distributions of 
(i) the total number of collisions and (ii) the ordering time T Q . We show that several 
features of these distributions are universal. 
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1. Introduction 

Exactly solvable models of interacting particle systems are important both in equilibrium 
and out-of-equilibrium statistical mechanics. Such solvable models, apart from their 
pedagogical virtues, also offer important insights and intuitions of the underlying 
more complex physical phenomena. Besides, they also allow to understand better the 
limitations of the approximate methods one generally uses in treating many-particles 
systems (e.g., the Boltzmann equation). 

One such very useful and instructive model was introduced a few decades ago [1, 2] 
and is usually referred to as the Jepsen gas [3]. It consists of identical hard-point particles 
that undergo binary elastic collisions on a line. At these instantaneous collisions the 
particles exchange their velocities, while in-between collisions they move freely. Due 
to the simplicity of the dynamics, this system admits analytical treatments for various 
imposed conditions, and therefore triggered quite a lot of interest in the past. The work 
of Jepsen [3] was followed by Lebowitz et al. [4, 5, 6], McKean [7], and Keyes [8] who 
refined and extended the calculations of the stochastic properties of a 'test' gas particle, 
including its asymptotic diffusive-like behavior and a comparison with the results of the 
Boltzmann equation approach. There was a recent surge of interest in this model in the 
context of the "adiabatic piston problem" [9, 10, 11, 12, 13], including the description 
of the stochastic dynamics of the piston-particle [14] and the analytic calculation of 
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the energy transfer and heat flux throughout an inhomogeneously-prepared, out-of- 
equilibrium configuration of the gas [15]. Later on, the Jarzynski theorem [16] was 
illustrated for the case of uniform expansion or compression of the gas [17]. Jepsen gas 
also turned out to be useful in the study of spin transport processes in the nonlinear 
a-model [18, 19]. 

In the context of biological evolution, a simplified version of Eigen's quasispecies 
model [20], namely, the shell model [21, 22, 23] corresponds to the free expansion of 
a Jepsen gas of N particles, where initially there is one particle at each position —k 
for k = 1,2, ...,N with a positive velocity drawn independently from a position 
dependent probability density function (pdf) 4>k{Uk)- A further simplification treating 
the velocities as independent and identically distributed (i.i.d.) random variables with 
a common k independent pdf (i.e., <pk = 0?V7c) leads to the i.i.d. shell model, for which 
several asymptotic properties can be computed exactly [24, 25] . In particular, the statis- 
tical properties of the piston-particle (the rightmost particle or the leading genotype in 
the biological language) exhibit universal properties. Its velocity distribution function 
at intermediate times (1 <C t <C iV 7 , where 7 is related to the tail of the velocity dis- 
tribution) has universal scaling behavior of only three varieties, depending exclusively 
on the tail of 4>{U) [25]. The associated scaling functions are different from the usual 
extreme-value distribution of uncorrelated random variables, and this difference is due 
to the dynamically-built, collision-generated correlations between the particles at finite 
times. These correlations are also responsible for the fact that the statistics of the pis- 
ton collisions is not Poissonian. Indeed, both the mean and the variance of the number 
of collisions the piston-particle undergoes increase logarithmically with time, but with 
different prefactors (that are also universal). 

Despite this rather rich history of the Jepsen gas, there are other interesting natural 
questions to which detailed answers are missing. Consider, for example, the fact that 
the dynamics of the freely expanding Jepsen gas provides a "natural" sorting algorithm 
for the velocities of the particles. It leads to an asymptotic "fan" state (see figure 1) in 
which these velocities are completely ordered from left to right in an increasing fashion. 
Two questions arise naturally from this ordering of the Jepsen g cLS £1X6 

A. What is the total number of collisions the particles undergo? 

B. How long does it take for the gas to reach the "fan" state? 

In this paper we provide analytical answers to these two questions. 

More precisely, we consider the free expansion of a Jepsen gas of N hard-point 
particles of equal mass on the infinite real line (—00,00). The initial positions of the 
particles are drawn independently from a common pdf ij)(X) and a random velocity 
drawn independently from a common pdf 4>(U) is assigned to each particle. At 
subsequent times (t > 0) each particle moves ballistically according to its assigned 
velocity, and upon contact between two particles they undergo elastic collision which 
merely interchanges their respective velocities. We assume both the set of positions and 
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Figure 1. (Color online), (a) A realization of free trajectories for N = 4 particles. 
The trajectories are labeled according to the order of their starting positions from 
left to right, Y\ < Y-x < I3 < Y4. The slopes {U\, U2, U3, U4} of the trajectories 
with respect to the "TIME" -axis represent velocities associated with them, (b) Actual 
particle-trajectories. After a ordering time the system reaches a "fan" state where the 
velocities are completely ordered from left to right in an increasing fashion and there 
is no further collision. The total number of binary collisions N c = 3 



the set of velocities are continuous variables (or at least one set), so that the collisions 
are always binary, and there can be at most one binary collision at one instant of time. 
Clearly, in each collision the velocities of the colliding particles get ordered such that 
after the collision the particle on the right acquires the larger of the two velocities and the 
one on the left gets the smaller of the two velocities. Therefore, after a certain ordering 
time the system reaches a "fan" state where the velocities of the particles are increasingly 
ordered from left to right. Once this "fan" state is reached, evidently, there cannot be 
any further collision. For any given initial realization of positions and velocities, the 
dynamics of future evolution of the gas is completely deterministic. Therefore, both the 
total number of collisions (denoted by N c ) and the ordering time (denoted by T Q ) are 
solely determined by the initial condition. Hence, N c and T Q are random variables in 
the sense that they differ from one realization of initial condition to another. In this 
paper we analytically compute their distributions. Our main results are: 

A. The probability Pn(N c ) of having N c collision, is completely independent of ip(X) 
and <f)(U) for all A ^ 2, and for large A it approaches a Gaussian form around its 
mean (N c ) = A(A - l)/4 with a variance (A c 2 ) - (A c ) 2 = N(N - 1)(2N + 5)/72. 

B. When the ordering time T Q is suitably scaled with A as r Q = T /[bN 2 ], — where b is 
a nonuniversal scale factor which depends explicitly on ipiY) and 4>{U) as given by 
(16), — the limiting pdf of r Q in the scaling limit A — > 00, T D — > 00 while keeping 
r Q fixed, becomes universal, i.e., completely independent of i>(X) and 4>(U), and is 
given by the well known Frechet form f(r ) = t q ~ 2 exp(— l/r G ) that arises in the 
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extreme value statistics. 

The paper is organized as follows. In section 2, we compute the statistics of the 
total number of collisions, and in section 3, we compute the limiting distribution of the 
ordering time. Section 4 contains some concluding remarks. Most of the calculational 
details are relegated to the Appendices A-C. 

2. Total number of collisions 

In order to express the total number of collisions N c in terms of the initial condition, we 
label the particles as i = 1, 2, . . . , N from left to right, i.e., i — 1 is the leftmost particle 
and i = N is the rightmost one (see figure 1). Let Y{ for each % = 1, 2, . . . , N denote the 
initial position of the z-th particle such that Y\ < Y 2 < • • • < Yjv- Note that the labeled, 
ordered coordinates {Yj}'s are no longer distributed independently according to ij)(Yi), 
rather their joint pdf is given by 

N N-l 

^jointO^ Y N ) = N\ J] ip(Yi) H 9(Y 3 - Y j+1 ), (1) 

i=l j=l 

where 9{x) is the Heaviside step function. On the other hand, the initial velocities 
associated with the labeled particles, which we denote by Ui for i — 1, 2, . . . , N are i.i.d. 
random variables drawn from the common pdf 4>{U). 

Now, for a given initial configuration, the system is fully and uniquely described at 
all subsequent times by the set of the free trajectories [Yk + Ukt \ k — 1, 2, . . . , N and t ^ 
0} (see figure 1(a)). Note that the k-th free trajectory should not be confused with the 
actual trajectory of the k-th particle. Indeed, each of the particles travels along such 
a free trajectory until it collides with another particle; in such a binary collision the 
particles interchange their trajectories (see figure 1(b)). In terms of the free trajectories, 
the total number of collisions N c is just the total number of intersections among the N 
free trajectories. Two free trajectories can, of course, intersect at most once. 

Consider first two extreme situations. If in the initial configuration the velocities 
are already in the increasing order, i.e., U\ < U2 < • ■ • < (which happens with 
probability 1/N\), then there cannot be any collision at subsequent times as the gas 
evolves. Thus, in this case, the system is always in the "fan" state and iV c = 0. On the 
other hand, if the initial velocities are in the decreasing order, i.e., U\ > U2 > ■ ■ ■ > t/jv 
(which also happens with probability 1/N\), then each pair of free trajectories intersects 
once and there are ( 2 ) of them. Therefore, this configuration yields the maximum 
number of binary collisions, which is N c = N(N — l)/2. For any other realizations of 
the velocities, iV c lies between and N(N — l)/2. The total number of collisions N c is a 
random variable, N c G [0, N(N — l)/2], which differs from one realization of the initial 
configuration of the particle velocities to another, and whose statistical properties we 
address below. 

Since, Yj < Y 3 - for i < j, two free trajectories (Yi + Uit) and (Yj + Ujt) intersect if and 
only if Ui > Uj for % < j. Therefore, for a given realization of the initial velocities {Ui}, 
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the total number of collisions (total number of intersections among N free trajectories) 
can be expressed as 

AT-l N 

^ = E E w-^)- ( 2 ) 

i=l j=i+l 

Obviously, N c is independent of the set of initial positions {Yj} of the particles and their 
distribution. Moreover, we will show below that Pn(N c ) is also completely independent 
of the velocity distribution <f)(U) and is solely determined by N. 

The probability P/v(A r c ) of having N c total number of collisions can be formally 
expressed as 



Pn{N c ) 



5 



N-l N 



]>(£/*) dC/* 



(3) 



i=l 



where S[n] with integer n, is discrete delta function: 
Let us consider the change of variables [26] 

(j){U)dU, i.e., dm = (j)(Ui) dUi for i 



U; 



1 and 5[n] = for any 

L,...,iV. (4) 



Obviously is a monotonically increasing function of C/j, and therefore 9(Ui — Uj) = 
9{ui — Uj). Moreover, since 4>(U) is normalized to unity, the variables Uj vary from to 
1, and so (3) becomes 

-i »i r N—i n 

P N (N C ) = d Vl ... dy N 5 N C -^Y1 6 ^ ~ u ^ ' ( 5 ) 

J o [ i=l i=i+l 

where the velocity distribution 0(C/) simply drops out, meaning that Pn(N c ) is universal, 
i.e., is the same as if the "new" velocities {ui} are drawn independently from an uniform 
distribution over [0, 1]. 

The mean number of total collisions (N c ) is straightforward to compute. Using the 
change of variables in (4), it is trivially checked that (9(Ui — Uj)) = 1/2. Therefore, 
taking the average in (2) yields 

N ~ l N N(N-U N 2 

TO = E E ~ U >))= 4 ~ — asiV^oo. (6) 

i=l ,3=1+1 

The calculation of the variance implies several steps that are indicated in Appendix 
A. One obtains finally: 



a 



,9, , , 9 N(N-1 27V + 5 N 3 
iV r 2 - (N c 2 = — i as TV 



OO. 



(7) 



As given by (2), N c is a sum of N(N — l)/2 random variables that are correlated. 
However, each term in (2) is correlated with only 2(iV — 2) other terms, since two 
different terms are correlated only when they have one velocity in common. Therefore, 
near the mean (N c ) and within a region 0(a), the probability distribution Pn(N c ) for 
large iV has a Gaussian form (see Appendix B), 



Pn(N c 



2txo 



exp 



2a 2 



(8) 



Number of collision and ordering time in the Jepsen gas 



7 




Figure 2. (Color online). Probability density of the total number of collisions (scaled 
by the mean and the variance) for N = 25, 50, and 100 (<>, V, and + respectively as 
symbols), as compared to the Gaussian p(x) = exp(— cc 2 /2)/v / 27r (drawn with ). 
We scaled the numerically-obtained data by the mean (N c ) and the variance a 2 that 
are given as functions of N by (6) and (7) respectively. It is difficult to notice any 
difference between the three data sets and the Gaussian. 



One can check numerically that this is actually already an extremely good approximation 
for N as small as 25 and for practically all the values of N c , as illustrated by the figure 2. 



3. Ordering time 

In order to compute the statistics of the ordering time T Q , it is convenient to relabel 
the particles according to the decreasing order of the velocities (see figure 3). Let 
{Vi, Vii ■ ■ ■ , Vn} denote the decreasingly ordered set of the initial velocities, i.e., 

Vi =m&x(Ui,U 2 ,...,U N ), 

V 2 =m8x({U 1: U 2 ,...,U N }\{V 1 }), 

V 3 =max({U lt U 2 ,...,U N }\{Vi,V 2 }), 



V N = min(L r i, U 2 , . . . , U N ), 
so that V\ > V 2 > ■ ■ ■ > V^. The joint pdf of these ordered velocities is given by 



0joint(Vl, V 2 ,..., V N ) 



N 



X 



N-l 



I W - v i+1 ] 



(9) 
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Figure 3. (Color online). The same realization of the free trajectories shown in 
figure 1. However, now the trajectories arc labeled according to the decreasing order of 
the velocities V\ > V2 > V3 > V4 and the initial positions {X\, X 2 , X 3 , X4} are indexed 
according to the order of the velocities. Note V\ = U2, V2 = U4, V3 = U\, V4 = U3 and 
X\ = Y 2 ,X 2 =Y 4 , X3 = Yi, X 4 = y 3 , while comparing with figure 1. 



Let Xi for each i = 1,2, ... ,N denote the initial position of the particle having the 
initial velocity VJ. Note that when we label the particles according to the order of their 
velocities, their initial positions {Xj}'s are no longer ordered on the line, but are i.i.d. 
random variables drawn from the pdf ip{X). 

As before, for a given initial condition, the system is fully and uniquely described at 
all subsequent times by the set of free trajectories {Xi + V(t \ i — 1, 2, . . . , N and t ^ 0}, 
where V± > V% > . . . > Vjv- The "fan" state is reached when all the free trajectories 
become completely ordered according to the velocities, i.e., X\ + V\t > X 2 + V 2 t > ■ ■ ■ > 
Xn + Vftt at all t > T a (see figure 3). 

Let F N (T) be the cumulative probability distribution of the ordering time, i.e., 
F N (T) = Prob[T Q < T]. It immediately follows that F N {T) = Prob[Xi + V X T > 
X 2 + V 2 T > ■ ■ • > X N + VnT], given that V\ > V 2 > ■ ■ ■ > Vn- This probability can be 
formally expressed as 

F N (T) = (T[0([X i - X i+1 ] + T[Vi - V i+ i]) \ , (10) 

where (. . .) denotes the averaging over both the i.i.d. initial positions that are drawn 
from the common pdf ip(X), and the initial ordered velocities, whose joint pdf is given 
by (9). Note that, -Fjv(O) f° r an Y finite iV is non-zero, 

IN-l \ 

mo) = ( n K Xi ■ Xm ) ) = m (n) 
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and is universal. This represents the probability that the initial positions of the particles 
are also ordered according to their velocities, i.e., X\ > X 2 > ■ ■ ■ > X^, so that there 
cannot be any collision. The pdf of the ordering time is simply 

5(T ) , dF N (T Q ) 



f N (T ) m , ^ , 

where the extra weight at T Q = vanishes in the limit N — ► oo. 

Let us introduce the variables x^ = X i+ i — Xi and Vi = Vi 
terms these variables we rewrite the expression (10) as 

IN-l 



Fn(T) 



where v,- > for i — 1, 2, 



n 



i-e(x t )e{^- 



(12) 



Vi + i for brevity. In 



(13) 



N-l. 



Since {Xj}'s are drawn independently from ip(X), two different variables Xi and Xj 
with i 7^ j are correlated iff \i — j\ = 1. Therefore, for large N, treating {xj}'s as i.i.d. 

/oo 
4-'(X) if){X + x) dX, is a rather good 
-oo 

approximation, which becomes exact in the limit N — >• oo. On the other hand, in the 
variables — Vi — V i+ i (i — 1, 2, . . . , N — 1), the ordered velocities {V^} are correlated, 
as can be seen from their joint pdf given by (9). Therefore, the random variables 
ti = 9(xi)[xi/vi] with i = 1, 2, . . . , N — 1, are not independent either. It turns out, 
however, that the limiting distribution of the ordering time T Q = max(t 1 ,t 2 , • • • , ijv-i)j 
when suitably scaled with N, has the well known Frechet form (see figure 4) that arises 
in the extreme value statistics of i.i.d. random variables drawn from a common parent 
distribution with a power-law tail. To show this, we formally expand the product in 
(13) as a sum of terms 



F N (T) 



N-l 

E 

i=i 

N-2 I* 

E 

i=l j=i+l 



Xi)9 



f 



+ (-1)" 



+ ■ 



N-n N-n+1 



E 




(14) 



We show in Appendix C that T scales as N 2 when both T and N are large. 
Moreover, as detailed in Appendix C, considering first n N, and then taking the 
scaling limit T — > oo and — > oo, but keeping T/N 2 fixed, the n-th term (counting 
from the "zero"-th term, which equals unity) in (14) becomes 

N-n N-n+1 N-l 



'2lE e ■■■ e <n 



h=l i2=h+l 



l+1 \f=l 



0{x lv )9 



N 2 t' 



1 

n! 



(15) 
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Figure 4. (Color online). The pdf of the ordering time T Q , plotted using the scaled 
variable r Q = T /\bN 2 ], where b is given by (16). The points (©) were obtained 
by numerical simulation with iV = 1000, averaging over 10 8 realizations of initial 
configurations where particles were distributed uniformly and independently in the 
interval [0,1], and the velocities were drawn independently from a Gaussian pdf 

<j>(V) = cxp(-y 2 / 2 )/V27r. Equation (16) gives b = 1/[12 Vtt]. The solid line ( ) 

represents the Frechct pdf /(r Q ) = t ~ 2 cxp(— 1/t g ). The inset displays the same pdf-s 
on a logarithmic scale. 



where 

b = 

with, recall, 

g{x) 

Thus, (14) yields 



x g(x) dx 



x 



\V) dV 



ij(X)^)(X + x) dX 



(16) 



(17) 



lim F N (bN 2 



E 

n=0 



-1] 



IV. 



exp 



The limiting pdf of the scaled ordering time r Q = T /[bN 2 ], — where the scale factor b is 
nonuniversal and is given explicitly by (16), — is therefore given by the universal function 
f{ T o) = To" 2 ex p(— l/ r o)- Figure 4 compares the results of the numerical simulation of 
the gas with this limiting pdf. 
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The Frechet cumulative distribution (18) usually emerges as the cumulative 
distribution of the maximum of a set of N 2 /2 i.i.d. random variables {xj} each drawn 
from the common pdf having a power-law tail p(r) ~ 2b/ V 2 . In the context of Jepsen 
gas, let us consider two free trajectories (Xi + V{t) and (Xj + Vjt) chosen at random. 
The ordering time between these two free trajectories has the pdf 



p(r k 



S 



Tij - e 



Xj — Xi 
V,. - V, 



Xj — Xi 
Vi - V, 



2b 



for large r»,-, (19) 



where b is given by (16). The ordering time T Q of the full system is clearly the maximum 
of the ordering times {t^} between each of the N(N — 1)/2 (« N 2 /2 for large N) pairs of 
free trajectories. Therefore, it turns out simply that in the N — > oo, the ordering times 
between different pairs of free trajectories become uncorrelated, and can be treated as 
i.i.d. random variables drawn from a common pdf having power-law tail p(r) ~ 2b j V 2 . 



4. Concluding remarks 

In this paper we have computed the statistics of the total number of collisions of a freely 
expanding gas of N hard-point particles, and found that the mean is 0(iV 2 ). On the 
other hand, the variance is 0(N 3 ), due to the correlations between different particle 
collisions, without which the variance would also have been 0(iV 2 ). However, despite 
these correlations, the probability distribution of the total number of collisions near the 
mean is very well approximated by a Gaussian form. 

In the context of biological evolution of quasispecies, the evolution time T e — which 
is defined as the time at which the rightmost particle undergoes the last collision, i.e., 
the free trajectory with the largest slope (with respect to the time-axis) becomes the 
rightmost trajectory — was studied recently [21, 27]. It was estimated that its pdf has 
a power-law tail p{T e ) ~ T~ 2 . The ordering time studied here represents obviously 
the upper bound to this evolution time, T e ^ T Q . We have found that the limiting 
distribution of the ordering time, when suitably scaled with N, becomes universal and is 
given by the Frechet form. This Frechet form usually appears as the limiting distribution 
of maximum of a set of i.i.d. random variables drawn from a common distribution with 
a power-law tail. Our result here provides another mechanism of generating Frechet 
form as a limiting distribution. 

Finally, while we have studied the statistics of the total number of collisions, it 
would be interesting to study the number of collisions of the particles N c (t) up to time 
t. Our results correspond to the limit t —>■ oo. In this context we point out that 
recently the collision statistics of a tagged particle in a rf-dimensional hard-sphere gas 
at equilibrium was investigated using Boltzmann equation approach [28] and was found 
to be non-Poissonian. It would be interesting to verify this conclusion for a tagged 
particle in the l-d Jepsen gas for which an exact result may be possible to obtain. 
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Appendix A. Variance of the total number of collisions 

Subtracting the mean given by (6) from (2), then taking the square and the average 
gives 

, N-l N N-l N 



a' 



N c - (N e ) 



i=l j=i+l k=l l=k+l 



where 



2 0(Ui - Uj) - 1 = ±1. 



(A.l) 



(A.2) 



Note that, ( 8ij ) = 0, (s?.) = 1 and (s ijSkl ) = A{6(U l - Uj) 6(U k - U t )) - 1. 

The total number of terms in the summation in (A.l) is [N(N — l)/2] 2 . These 
can be grouped according to the correlation functions (siys^)'s that are of similar kind. 
This can be conveniently represented by diagrams as shown in Fig. Al. Then using the 
change of variables (4), it is easy to compute the correlations (6(Ui — Uj)6(Uk — U)) and 
hence (s^Sfez), corresponding to each of the diagrams (a)-(f) in figure Al. 
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Figure Al. Diagrammatic representations of the different types of terms appearing 
in (A.l) for the computation of (N%). In all the diagrams j > i and I > k. The solid 
lines connect two indices having equal value. One has: (a) i = k and j = I; (b) i = k 
and j I; (c) i ^ k and j = I; (d) i = I; (e) j = k; (f) i ^ k and j ^ I. 



Diagram (a): (sy-s w ) = (s 4 2 ,) = 1. 
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The number of such terms in (A.l) is equal to number of ways of choosing two 
indices out of N, which is 

'N\ N(N - 1) 



2 J 2 
Diagram (b): (%s w ) = (sijsn) = 1/3. 

The number of such terms in (A.l) is equal to number of ways of choosing three 
indices {i,j, 1} out of N with allowing permutation between two of them {j, I}. 
This is given by 

2x /JV\ _ N(N-l)(N-2) ^ 

Diagram (c): Same as for diagram (b). 

Diagram (d): (%s w ) = (sySfci) = -1/3. 

The number of such terms in (A.l) is just the number of ways of choosing three 
indices {k, i, 1} out of N, which is 

'N\ N(N -1)(N -2) 



3 ) 6 
Diagram (e): Same as for diagram (d). 

Diagram (f): (sijSki) = (sij)(ski) = 0. These terms do not contribute to the sum 
in (A.l). 



Therefore, finally (A.l) yields 
1 
1 



n ) + 2. 1 -. r 

2 3 V 3 



AT(AT-l)(2Ar + 5)_ 



72 

Appendix B. Probability distribution of the total number of collisions 

Let us consider the deviation of the total number of collisions given by (2) from it mean 
given by (6): 

1 JV-l N JV-l N 

m = n c -(n c ) = -J2Y, s ^ - u i) = EE s * (B-i) 

i=l j=i+l i=l j=i+l 

where sgn(x) = 26 (x) — 1, and Sij = (1/2) sgn(£/j — Uj). Obviously (M) = and we 
have already calculated the second moment in Appendix A, which is 

(J|05g , = *(*-l)(iW + 8) w g 

\ / 72 36 & K ' 

Assuming the pdf of the velocities to be symmetric, i.e. <f)(—U) = <p(U), and noting 
that sgn(— x) = — sgn(x), it can be easily shown that the probability distribution of 
M is symmetric about zero. On the other hand, we have shown in section 2, that the 
probability distribution of N c , and hence that of M, is independent of the velocity 
distribution. Therefore, the distribution of M must be symmetric for any velocity 
distribution. This implies the vanishing of all the odd moments, i.e., 

(M 2n+1 ) = 0, for n = 0,1,2,.... (B.3) 
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By the same argument, one can also show that average over any product of odd number 
of Sf.^'s is zero. (M 2n+1 ) is just a sum of such terms. 
The even moments are given by 

JV-l N N-l N I 2n \ 

<m 2 "> = £ £ •••£ £ (n s «»). < B - 4 ) 

for n = 1, 2, 3, . . .. We first break the average as product of averages over n different 
factors of pair of Si i's, and then sum over the indices independently under each average 
of factors, each of which gives the second moment (M 2 ) = a 2 making it ~ a 2n . This is 
0(iV 3n ). The corrections to this come from the terms where the sets of indices between 
different factors are not independent. If there is one common index between two different 
factors, we will have to take the average of them together. It will reduces one summation 
and therefore is 0(iV 3ri " 1 ). The number of ways of grouping 2n objects into n pairs is 
clearly (2n)\/[n\ 2 n ]. Therefore, 



(M 2n ) 



(2n)! 
n\2 r 



a 



2n 



+ O (iV 3 "" 1 ) , for n 



2,3, 



Let us now consider the scaled variable 

r M 
z = hm — . 

n^oo a 



Since a 2n is 0(A r3n ), we have 

(2n)! 



n\2 r 



and (z 2n+1 ) = 0, for n 



0,1,2,3, 



The characteristic function of the probability density of z is then 



= £ 



n=0 

oo 

£ 

n=0 



(2n)\ 

:-i) ra 

n\ 



X 2n ^ 



exp 



2 



The inversion of this Fourier transform gives the pdf of z as 



p(z) 



which says that the scaled variable [N c — (N c )]/a has a Gaussian distribution. 



(B.5) 



(B.6) 



(B.7) 



(B.8) 



(B.9) 



Appendix C. Evaluation of the terms in the series expansion of F^(T) 
First term: 

Let us compute the average (9(xi) 9(xi/T — Uj)), in the first term (counting from the 
"zero"-th term, which equals unity) in (14). The pdf of Vi is given by 

N\ 

Pi(v) = 9(v) 



(JV- i-1)! (i-1)! 



x / dV(j)(v) i-$(v) (piy + v) $(v + v) 



N-i-l 



i-1 



(C.l) 



Number of collision and ordering time in the Jepsen gas 



15 



in which 



<f>(w) 



<f>(V) dV, 



i.e. 



-<j>(V) dV. 



(C.2) 



The integrand in (C.I) merely specifies that (N — % — 1) velocities are smaller than V, 
(i — 1) velocities are larger than (V + v), two velocities are V and (V + v) respectively, 
such that their difference is v. The combinatorial prefactor gives the number of such 
arrangements, and we finally integrate the integrand over all possible values of V. It 

POO 

is checked that this pdf is normalized, / pi(v)dv = 1. Using (C.I) we first take the 

Jo 

average over Vi, which for alH = 1, 2, . . . , N — 1 gives 



e(xi 



e(xi) m 



(N 



1 ! i-1 ! 



x 



dV4>{V) 



N-i-1 



1 - $(V) 

®(V + Xi /T) 



When T > iV, 



,1, v ; 

T 



W + J 



(C.3) 



(C.4) 



Note that the pdf of Xi as given by (17) is independent of the index %. Substituting 
(C.4) in (C.3), then taking average over Xj, and finally summing over % yields 

g (ewe( Ei_„,)). 6 fci) + o(^, (c,) 



i=l 



T 



with b given by the (16). 

Equation (C.5) suggests that T scales as iV 2 , which indicates that the natural 
scaling limit corresponds to T — ► oo and N — > oo, while keeping T/N 2 fixed. In this 
limit the term (C.5) becomes 

N-l , 



1=1 



Second term: 

Let us now consider the second term in (14) involving summation over i and j. We 
substitute T = N 2 t', and rewrite it as 

N-3 N-l 



i=l J=j+2 
W-2 



X j 



N 2 t' 



vA6{ Xj 



N 2 t' 



— f i I P ( TTT^ 7 - f i+1 



1=1 



N 2 T 



N 2 t' 



(C.7) 



(C. 



Now in the first term (C.7), the variables Xi and Xj are uncorrelated, and each of them 
has the same pdf (17). Following the similar reasoning that we used to write down the 
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pdf of the single variable Vi in (C.l), we can also write down the joint pdf of Vi and Vj 
for j ^ i + 2 as 

PuM = {N _ j _ 1)l{j _ i _ 2)l{i _ 1)l 

N-j-l 



(N-j-l)\(j-i-2)\(i-l)\ 

poo . 

x / dV 4>(V) / d7>(7') 1 - 



0(V - V - uj) <f>{V + «<) 0(V + Vj ) 

i-i-2 r 



(C.9) 



where $(u>) is given by (C.2). It is checked that the joint pdf is normalized to unity, 

poo poo 

i.e., / dvi / dvj pij(vi,Vj) = 1. Now using (C.9) we first compute the average over 
Jo Jo 

Vi and Vj in (C.7), then we expand the result in Taylor series (assuming large N), and 
then average over Xi and Xj. Finally, we sum over i and j, and take the limit iV — > oo. 
Noting that 

1 2 

, (CIO) 



/CO /'OO 1 r POO 

dV<f> 2 (V) dV'cf> 2 (V') 9{V - V) = - / <f> 2 (V)dV 
-oo J — oo J —oo 



we find 



JV-3 JV-l 



J™ E E (o&mxi) 9 



1 

2! 



(en; 



i=l j=i+2 

where 6 is given by (16). The term (C.8) goes to zero in the limit N — > oo. 



JTie n-t/i ierm: 

Following the exactly same steps, we can evaluate the n-th term of the expansion of 
Fn(T) in the expression (14). We first assume that n <C N, and later take the limit 
N — > oo. We again write the sum as 

N-2n+l N-2n+3 N-l 



E E ••• E n 

il=l «2=n+2 i n =i„_i+2 \i/=l 



N 2 t' 



Vi 



+ ^remaining termsj , 



(C.12) 



(C.13) 



such that (C.12) contains only the terms with z& ^ i-k-x + 2 for all k = 2, . . . , n. As 
before, the variables {x iv , v — 1,2, ... ,n} in (C.12) are uncorrelated, and the joint pdf 
of {v iv } can be written as 



Ph,i2,...,i„( v in v i2i • • • ■> v i n ) 



n«<< 

y=X 

oo 



C N (i 1 ,i 2 , ...,i n ) 



x / dV (1) 0(V (1) ) / dV (2) 0(V (2) ) 



dV( n) 0(y (n) ) 



'n-l 



\{o(v {k) -v {k+l) ) 



.k=X 



oo 

JV-in-1 
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U°( v ^)- v w-0 

k=2 

n 



l<t>(y (k) +v ik ) 



,k=2 



k=l 

ifc— ife-1— 2 



ii-1 



(C.14) 



where $(u>) is given by (C.2) and where the combinatorial prefactor Cn gives the number 
of ways of arranging N velocities as in the above integrand, 

C N (h, z 2 , . . . , z„) = N\ - 1)! - z fc _i - 2)!] (N-i n - 1)! J . (C.15) 

It is again checked that the joint pdf (C.14) is normalized, i.e., 



l 2 , 



^w^,^,...,^) JJd^ = 1. 



(C.16) 



Using (C.14) we first take the average over in (C.12), then we expand the result 
in Taylor series (assuming large N), and then average over {x^}. Finally, summing over 
all the indices {i u }, and taking the limit iV — > oo, (C.12) yields 

1 



(r') n 



x g(x) dx 



o 



ldV [k] <p 2 (V {k) ) 



,k=l 



n-1 



[[e(v {k) -v {k+1) ) 



k=l 



(C17) 



(C.18) 



Now, the term inside the first square brackets of the integrand in (C.18) remains 
unchanged under permutations of the variables {V(fc)}. While the term inside the second 
square brackets of the integrand changes under permutations, summing over all possible 
permutations yields unity. However, note that {V(fc)}'s are just dummy variables in the 
multiple integral (C.18), so it must remain unchanged under any of the n\ permutations 
of these variables. Therefore, the integral (C.18) equals 



1 



n 2 (vfo) 



k=l 



1 



cj) 2 (V) dV 



(C.19) 



Therefore, in the limit iV — > oo the expression (C.12) becomes [see (C.17), (C.18)], 

7V-2n+l N-2n+3 N-1 



>- E E - e <n 



0M9 



• ] j. 



N 2 t' 



Vi 



, (C.20) 



ii=l t2=h+2 i n =i„~i+2 \u=l 

where b is given by (16). The remaining term (C.13) goes to zero in the limit N — > oo. 
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